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1.  Introduction 


ed  °?e  Fas,'Fiew  program  (ffp>  «*>  - 

.ions,  based  on  geonrertcal^lf 

ground  impedance,  and  diffraction  because  of  refracb ’Ptl°n’  r'fraC,‘°n'  acoust'oaily  complex 
version  uses  the  output  of  SCAFFIP  to  build  mat  r  °Ver  fla(  terra,n-  The  Microsoft  Windows 
range.  The  attenuation  matrix  can  be  use  >s  Z  B „TT  «**"“»  a"d 

such  as  CASTFORM,  JANUS,  and  BASIS.  P  '  'C  OT  acoustic  sens»r  simulations 

s“opiX„”i!!0„s  Jur “eC«rV,3,Tdi°n 0ngrIIy  deVe‘°Ped for “nde™a«or 

gation  in  the  atmosphere.  The  FFP  developed  bv  R  /  ,L?  “  “  '  |d|  adapted  ^  FFP  <°  Propa- 
tion.  If  each  layer  in  the  atmosphere  is  viewed  as  an  o  t  “  propa8ati°.>  matrix  formula- 

atmosphere  can  be  constructed.  Multiplying  each  mar  T  *  m!l'nX  ^  ^ ,ayer  in  the 

represents  how  an  acoustic  signal  will  he  f r  .a  ■  *  together  results  in  a  new  matrix  that 

a  Bessel  Function  Transform  of  the  problem  is  tak  **  "  pa°pa8ates  throu8h  the  atmosphere.  Next, 

model  on  a  PC-based  computer  platform  runnLgMcrTsoftZdr  t" TV 

cases  (m  Appendix  A)  show  the  structure  of  the  input  files  used  by  WSCAFFtP  a  Sa”P'e  ** 
a  basis  to  verify  their  installation.  The  software  0  ■  ?  ^  WSCAFFIP  and  give  the  user 

mail.  The  current  version  of  WSCAFFIP  dnP<  ?  ♦  § “  SVaiable  on  several  media,  including  e- 

nmst  copy  the  files  to  his  or Z  feZZ rZ  I'T  “ pr0gram’  S0  **  -r 

which  assists  in  copying  the  files  to  the  correct' 

procedure  for  installing  WSCAFFIP.  Sectl0n  4- 1  Slves  the  step-by-step 


2.  Propagation  Theory 


2.1  Speed  of  Sound 

Meteorological  phenomena  can  have  a  significant  ^ 

meteorological  variables  affecting  the  speed  of  sound  °"  '  re°e'Ved  S°U"d  field  S°me  of  the 

velocity,  and  humidity.  To  observe  the  effect  nf  ^  are  pressure.  temperature,  wind 

observe  the  effect  of  each  meteorological  variable,  consider  each  one 


independently  and  examine  the  equation  for  the  speed  of  sound  in  air.  The  value  c,  according  to 
Laplace’s  adiabatic  assumption  for  air  as  an  ideal  gas,  [5]  is 


(2.1) 


in  which  y  is  the  ratio  of  specific  heats,  R  is  the  universal  gas  constant  equal  to  8314.16  J/(kg  K), 
and  Mis  the  molecular  weight  of  air. 


The  presence  of  water  molecules  alters  the  sound  speed  by  lowering  g  and  decreasing  M.  The 
decreased  M  dominates  so  that  the  overall  effect  of  increasing  humidity  is  an  increasing  sound 
speed.  These  changes  can  be  quantified  as 


and 


7  +  h 
5  +  h 


(2.2) 


M=29-  11  h 


(2.3) 


in  which  h  is  the  fraction  of  water  molecules  in  air.  As  the  amount  of  water  in  the  atmosphere 
increases,  the  molecular  weight  of  a  parcel  of  air  will  decrease  since  the  molecular  weight  of  a 
water  molecule  is  less  than  diatomic  nitrogen  and  diatomic  oxygen.  As  a  result,  the  sound  speed 
increases  as  the  fraction  of  water  molecules  in  the  air  increases.  In  order  to  calculate  the  fraction 
of  water  molecules  in  air,  the  Goff-Gratch  equation  must  be  used  to  first  calculate  the  partial 
pressure  of  saturated  water  vapor,  Psat  at  temperature  T.  This  equation  is 


fp  ) 

(  t  A 

r  T  ^ 

x  sat 

=  10.79586 

1± 

i01 

±  5.02808  log10 

— 

lPo> 

It  J 

vToi> 

+  1.5047xl0-,(l-10-o,“l<T,T">-'1)  (2.4) 

+  0.42873x1  Cr'(lO‘-™”|MI-  Tl1 -l)- 2.2 195983 

in  which  7oi  =  273.16  K  and p0  =  1  atm  or  the  reference  pressure. 


After  the  value  for  Psa,  is  determined,  the  fraction  of  water  molecules  in  air  can  be  calculated  by 
the  following  relationship: 

jo-awif, 

P  V 

in  which  RH  is  the  relative  humidity  in  percent  and  P  is  the  pressure  in  atmospheres. 


The  magnitude  of  the  dependence  of  the  sound  speed  on  humidity  is  not  obvious.  To  understand 
the  degree  of  the  effect  of  humidity  on  sound  speed,  consider  a  particular  case.  At  20  °C,  the 
difference  in  sound  speed  between  0%  and  100%  humidity  is  2  m/s.  A  fluctuation  in  the  humidity 
of  this  amount  is  very  unlikely.  If  the  variation  in  humidity  is  reduced  to  a  change  of  50%  to 
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1 00%,  the  change  in  the  sound  speed  is  nnlv  1  m/c  -ru  r 

caused  by  changes  in  humidity  should  always  be  much  kss  thanl^ Ge0"  Sp6ed 

tuations  can  be  ignored,  ^enerahy,  humidity  flue- 

isctlculatedby6  W‘”d  Speed  °"  the  Speed  °f  S0Und  iS  a  Vec,0r  rela,i0"'  The  effective  sound  speed 

c«ff-c(T)  +  ux  «*(<,„- p-gj  (2  6) 

tzi  ;ne:he  abs:r  °f  wind  at  ^  1™^  of 

front  which  the  wind  bWs  C  9  ^  “  T*™  fiwn  ““  ■■  <»e  direction 

directions  are  relative o  LiT  ‘  "  ^  ™d  *  bl°™g  (~  Figure  I).  AU 


a“  "edt^^r  ^  S0“"d  ,S  3  “  °f  temperature 
through  the  atmosphere  The  degree  of  refract’  th  a°°UStlC  WaVe  t0  be  refracted  as  it  propagates 
speed  gradient  present  in  the a,™  her  f 2"  7*  Underg°eS  *  re'ated  tha 

wave  wiii  be  retracted  downw  T  f  th  a  l  “ **“  With  hei8ht. «*  '■““Stic 

will  be  refracted  upwtd  ****  d'CreaSas  w'tb  height,  the  acoustic  wave 


2.2  Fast-Field  Program  (FFP) 

The  propagation  of  sound  from  a  point  source  a,  the  origin  Js  given  by  the  classical 

-  j_£ 

c  dt 


v  P  2  ^  f --4^(x,y,z) 


wave  equation 
(2.7) 
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in  which  <5  represents  a  delta  function  source  of  unit  strength.  For  simple  harmonic  motion,  we 
may  replace  d/dt  by  -ico,  so  Equation  (2.7)  becomes  the  Helmholtz  equation 

V2p  -  k2p  =  -4 pS  (x,y,z)  (2.8) 

in  which  k  is  the  wave  number,  a j/c,  in  which  c  is  the  sound  speed  and  co  is  the  angular  frequency. 
For  the  FFP,  c  is  restricted  to  vary  only  in  the  z-direction. 

Transforming  Equation  (2.8)  into  cylindrical  coordinates  and  assuming  azimuthal  symmetry,  the 
Helmholtz  equation  becomes 


d2p  1  dp  d  f  1  dp)  2  2 

+  - ~T  +  k  p  =  -S(r)S(z-zs) 

or  r  or  oz\p  ozJ  r 


in  which  the  source  is  located  at  r  =  0  and  z  =  zs  and  p  is  the  density  of  the  medium. 


The  atmosphere  is  viewed  as  a  series  of  constant  sound  speed  layers  for  the  FFP  (as  shown  in 
Figure  2).  The  layers  in  the  atmosphere  are  bound  on  top  and  bottom  by  complex  impedance 
surfaces.  The  top  boundary  is  typically  modeled  as  an  infinite  half-space  with  constant  parameters. 
At  the  bottom  boundary,  the  atmospheric  layer  adjoins  a  partially  absorbing  surface  within  can  be 
represented  by  the  complex  acoustical  impedance  of  the  ground. 


Figure  2.  Graphic  representation  of  WSCAFFIP’s  view  of  the  atmosphere. 


To  reduce  the  dimensionality  of  Equation  (2.9),  a  zero-order  Hankel  transform  is  applied  with 
respect  to  the  range  variable  r.  We  define  the  transform  pair: 


P(K,z)=  f0  P(r > z)J 0(Kr)rdr 
p(f,  z)  =  fQ  p(r,  z)J0  ( KT)KdK 
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Applying  the  first  transform  to  Equation  (2.9)  results  in 


Jp  +  \k2(z)-K2\p  =  -2S(z-zJ  . 

Since  p„  du2  Idt  =  dp/dz ,  we  now  have  the  following  coupled  equations: 


(2.11) 


du2  _  i  d2p 

dz  copo  dz2 


~  _  i  dp 
cop0  dz 


(2.12) 


(2.13) 


in  which  ft,  is  the  mean  air  density  and  S,  is  the  transformed  particle  velocity  in  the  2-direction. 
The  delta  function  is  the  source  tenn  producing  a  discontinuity  in  u.  at  a  height  of  2,.  In  order  to 

thtLm  rc!hUl,i"8  T‘i0"S’ Ue  6<  a1'  Hl  USed  “  ana'0gy '°  3  transmission  line  that  results  from 
the  form  of  the  transfonned  equation.  The  transformed  equations  have  a  vety  similar  form  to  the 

e  egrapher  s  Equation  of  electrical  transmission  line  theory.  |4]  Looking  at  the  analogous 
telegrapher’s  equations,  6 


dl(z ) 
dz 


dV 

—  =  -ZI(z) 
dz 


-Y(z)V(z)  +  IoS(z-zs) 


(2.14) 


(2.15) 


there  is  a  similar  form  of  the  equations  if  the  shunt  admittance  Y(z)  is  made  a  function  of  a  and  the 
senes  impedance  Z a  constant,  with  a  current  source  at  2  >  2,.  The  equivalent  transmission  line 
configuration  to  Figure  2  is  shown  in  Figure  3. 
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Figure  3.  Transmission  line  analogy  to  the  atmosphere. 

Using  the  Lee  et  al.  analogy,  we  can  arrange  the  acoustic  problem  so  that  a  solution  can  be  calcu¬ 
lated.  We  make  the  analogy  by  representing  each  layer  in  the  atmosphere  by  an  element  in  a 
transmission  line  with  a  certain  characteristic  admittance  and  attenuation  constant.  The  admittance 
is  defined  as  the  inverse  of  the  impedance  or  the  inverse  of  the  sum  of  the  resistance  plus  the 
reactance  of  the  electrical  element.  The  admittance  of  the  element  causes  the  voltage  running 
through  it  to  be  attenuated  and  causes  a  shift  in  phase  of  the  signal.  Viewing  it  from  the  perspec¬ 
tive  of  the  acoustic  wave  problem,  the  acoustic  wave  undergoes  losses  and  phase  shifting  because 
of  refraction  and  spreading  of  the  acoustic  wave  as  it  propagates  through  the  atmosphere.  The 
equivalence  between  the  transmission  line  and  the  acoustic  wave  problem  is  illustrated  in  Table  1 . 
This  equivalence  can  be  performed  for  each  layer  of  the  atmosphere,  thus  constructing  a  transmis¬ 
sion  line.  The  top  and  bottom  boundaries  in  the  atmosphere  become  loading  admittance  elements 
on  each  end  of  the  transmission  line.  Now  the  problem  has  been  simplified  to  one  of  determining 
the  voltage  in  the  transmission  line  at  a  point  in  the  line.  This  is  a  well-known  process  in  electrical 
engineering. 


Table  1 .  Transmission  line  analog  relationships  to  acoustic  fields 


Acoustic  Wave 

Transmission  Line 

Fields 

Pressure  p(k3z) 

Voltage  V(z) 

Velocity  (o  (k,z) 

Current  I(z) 

Attenuation  Constant 

r = -  (&)  / 

r=  VzF 

Characteristic  Admittance 

Yc  =  -r/wp 

yc=  VFTz 
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The  zero-order  Bessel  function  in  Equation  (2  im  u 

Equation  (2.10)  can  be  expanded  via  Hankel  functions: 

Jo  i*r)  =  ~  [//W  (Kr)  +  HV)  (Kr  J 

T,  2  J  (2.16) 

The6  FFP  W  rustic  waves,  respective!, 

function  to  be  suppressed  with  the  additional1  argument  th^the  '  a**°WS  Por  tbe  ^rst  ^an^el  ^ 
contribute  significantly  to  the  final  result  The  1  T!  Inco™ng  acoustic  wave  will  not 
tion  for  large  arguments  is  as  follows:  ymptotlc  expansion  of  the  second  Hankel  func- 


2  e-u<^/4) 


"V  71 IK 


..... 

\!r 


Kr  »J 


(2.17) 


inverse  transform^^  fr°m  the  inte§rand  oft the 

(2.17)  into  the  inverse  transform  of  Equation  (2  1 01  ***  u  <  k°'  Substltutm§  Equation 

acoustic  pressure  equation  can  be  written  as  '  lng  the  far'field  approximation,  the 


(2.18) 


cal  integral  over  discrete  values  of^A^^  by  3  nUmeri' 


in  which 


^  r  «  --  0 


(2.19) 


—  f-  = 

to  NAk  ' 


(2.20) 

number  of  tenns.  ’  g  summation  to  be  terminated  at  a  finite 

A  problem  in  this  representation  of  the  solution  is  that  m  ) ;«  „•  , 

real  axis.  These  points  represent  •  h  *  ”  smgular  at  some  points  along  the 

the  current  FFP  uses  “extra  loss”  in  the  c^cuteT'08  Th^  3ttenU3tion-  To  avoid  this  difficulty, 
extra  loss  is  to  move  the  numerica|  ;  '  ”“he™tical  interpretation  of  using  this 

loss  from  the  solution  i„  an  approximate  ma„lr!  ^  We  rem°Ve  the  efftct  °f  «*» 

term  exp(a-),  in  which  ais  the  extra  loss  attenua(y  P  y'”S  con,Pu,ed  Pressure  by  the 

Proper  choice  of  this  “artificial  attenuation”  is  ess  t  iT  Np/m'  "  musl  be  nofecl  that  the 
from  the  code.  Cn  13  lf  meaningful  results  are  to  be  obtained 
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Another  problem  with  the  numerical  integration  is  the  number  of  points,  N,  used  in  the  summation 
of  Equation  (2. 1 9).  A  lower  bound  on  the  number  of  points  required  in  the  summation  is  |6] 


1C  V 

AT  __  ^  max'  max 
J  V  min  _ 

Tt 


(2.21) 


However,  in  most  cases,  this  lower  bound  is  too  large  to  perform  with  a  single  FFT.  The  problem 
is  how  to  have  enough  points  to  sufficiently  sample  the  wave  number  space  in  the  numerical 
integration.  If  Equation  (2.19)  is  rewritten  in  the  form 

S  =  (2.22) 

n=0 


this  summation  can  be  rewritten  so  that  the  single  summation  is  rearranged  to  a  double  summation 
of  the  form; 


P  N’-\ 

o  V  ^-'l2*(/>-l)m/pj  y  r  -4 2nnml{pN'\  (  , 

°  >  \L.L5) 

b~  0  «=0 

in  which  N'  =  N/p  and  p  is  an  integer  larger  than  1 .  p  is  the  number  of  panels  into  which  the 
original  FFT  has  been  divided.  Each  panel  contains  Appoints.  This  technique  allows  the  FFT  to  be 
performed  to  calculate  the  acoustic  pressure  with  range  at  the  height  of  the  receiver. 


2.3  Absorption  of  Sound  by  the  Atmosphere 

Viscosity,  heat  conduction,  and  molecular  exchanges  of  energy  basically  cause  losses  in  the 
medium.  In  the  1 9th  century,  the  mechanisms  of  viscosity  and  heat  conduction  were  the  only  ones 
suspected  of  causing  dissipation  of  sound.  Therefore,  they  are  presently  referred  to  as  classical 
absorption. 

In  classical  absorption,  if  one  represents  the  effect  of  absorption  by  a  factor  e~ar  in  which  r  is  the 
distance  of  propagation,  then  the  attenuation  coefficient  ac/  attributable  to  viscosity  and  heat 
conduction  is  given  by  Equation  (2.17)  from  Physical  Acoustics  XVII  [7] : 


a., 


5.57  xl(T9 


t/t„  r_ 

T  +  n0.4P/Po  • 


(2.24) 


The  units  of  ac\  in  Equation  (2.24)  are  nepers/meter,  in  which  Pa  is  the  reference  pressure  of 
1 .013  25  x  1 05  N/m2  (1  atm),  P  is  pressure  in  N/m2,  Ta  is  the  reference  temperature  of  293. 1 5K,  T 
is  the  temperature  in  Kelvin,  and /is  the  frequency  in  hertz. 

In  molecular  absorption,  energy  exchanges  at  the  molecular  level  include  rotational  and  vibra¬ 
tional  modes.  Analysis  of  the  rotational  mode  shows  that  its  representative  attenuation  coefficient 
is  proportional  to  ach  the  classical  attenuation  coefficient: 


^=4.16e-1687~V3 
a  a 


(2.25) 


when  29K  <  T  <  690K. 
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ca|ra!idli!oircuL^b^^tiori^areaddkiveeAdenl0rrrHte^  ^  ^  '°SSeS  a,tribu,able «»  classi- 
written  ^  dd’"Ve' A  Slmpllfed  empi™al  «>"»  of  the  equation  ca„  be 


«„  =  1.83x10 


-n  JrTr.P 


PIP 


(2.26) 


which  is  correct  within  2%  for  2 13K  <  T  <  373K 

a'SM  bC  °Tidered'  Because  ,he  -“■*«  -  cona¬ 
tion  absorption  is  given  by  ?  ” “  C°n,r,bUte  “  at« »”  coefficient.  The  vibra- 


av,hJ 


35c 


(2.27) 


(OSOMs’for oxy^Jd  08780<84"fro8Tn’ ^  ^  ^  m0k  fn,'‘io”  °f air  component  consid- 
ature  (2239. 1 K  for  oxygen  and  3352  OK  f  \  * the  Characteristic  vibrational  temper- 

Tin  m/s.  The/,  are  “7  C  "  **  °f  ^  *  *“« 

given  ^.Pe„dS  °n  tbe  nefative  humidity  and  .bnospheri!:  preLure.^h^daxadon  frequenci^m-e 


fr.o  =  —[24  +  4.04x10*  h—  + 

Poy  0.391+  h' 

fr,N  =  jr^(9  +  2%0he~Am{T°np'-^ 


(2.28) 


(2.29) 


\  - , 

The  total  attenuation  coefficient  is  then  ttip  cum  ^  ™  a 

attenuation  coefficient  of  T  =  20  °C  and  Rh  -  20°/  d  aT*  4  ^  3  log'Iog  plot  of  total 

contributions  to  attenuate  from  classmal  L r  T  ^  ”*«*'!*  11  Sh°WS  the 

oxygen,  and  the  total  attenuation  coefficient  attributa^  rdaXatl0n  ofnitro8en  and 

mechanisms.  coefficient  attributable  to  the  sum  of  these  three  attenuation 
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Figure  4.  Log-log  plot  of  sound  absorption  coefficient  versus  frequency  for 
sound  in  air  at  20  °C  at  1  atm  pressure  and  with  a  water 
vapor  fraction  h  of  4.676  x  10"3  (. Rh  =  20%). 


The  attenuation  coefficient  a  is  roughly  proportional  to  the  square  of  the  frequency.  Thus,  when 
frequency  doubles,  attenuation  will  quadruple.  The  attenuation  of  the  sound  wave  attributable  to 
molecular  and  vibrational  absorption  is  very  important  for  frequencies  greater  than  or  above 
250  Hz.  For  frequencies  below  250  Hz,  this  attenuation  does  not  contribute  much  to  the  total 
attenuation  of  the  sound  wave. 


2.4  Complex  Ground  Impedance 

Several  models  are  available  for  calculating  the  complex  ground  impedance.  The  impedance 
model  used  in  WSCAFFIP  is  the  Four-Parameter  Model  of  Attenborough  [8].  In  this  impedance 
model,  the  complex  normalized  characteristic  impedance  of  the  ground  is  calculated  by 

~V  .SW 

3Q  cop0  , 23 


in  which  q 2  =  W-"',  Sf  is  the  pore  shape  factor  ratio,  Q  is  the  porosity  of  the  ground,  <7  is  the  flo  w 
resistivity  of  the  ground  (inks)  rayls,  CO  is  the  angul  ar  frequency  (27$/),  p()  is  the  density  of  air  (1.2 
kg/m3),  and  kh  is  the  normalized  wave  number.  The  normalized  wave  number  is  computed  from 


Lzl 

Y 


\ 


Q  cop0 


(2.31) 


in  which  y  is  the  ratio  of  specific  heats,  Equation  (2.2),  and  Npr  is  the  Prandtl  number  (0.724). 
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e  parameters  Sfi  Q,  p,  and  n'  are  normally  varied  until  agreement  is  reached  between  impedance 
measurements  and  the  impedance  model  is  achieved.  However,  this  method  of  determining  the 
four  parameters  cannot  always  be  completed  if  time  or  resources  is  lacking.  To  aid  the  user  when 
the  four  parameters  are  unknown,  there  is  a  series  of  general  descriptions  of  a  variety  of  ground 

surfaces  from  which  to  choose.  The  user  has  the  option  to  enter  the  values  for  the  four  param¬ 
eters,  if  known.  F 

Some  ground  surfaces  have  a  layered  structure  that  results  from  the  gradual  deposition  of  material 
over  a  soil  base.  Some  examples  are  a  layer  of  snow  over  frozen  ground,  decomposition  of  or¬ 
ganic  material  over  day  or  sandy  soil,  or  a  well-plowed  pasture  over  clay  or  harder  packed  soil 
An  effective  impedance  Z(d)  can  be  calculated  for  a  semi-infinite  layer  of  impedance  Z2  covered 
by  a  layer,  depth  d,  or  another  material  of  impedance  Z, .  The  effective  impedance  is  given  by 


Z(d)  = 


z2  -  /Z,  lm(khd) 
[Z,-  iZ2  tan (kbd) 


Z, 


(2.32) 

in  which  k„  is  the  bulk  propagation  constant  in  the  top  layer  of  the  ground.  The  parameters  Z,  Z, 

and  k„  are  calculated  by  Equations  (2.30)  and  (2.31)  from  Attenborough's  Four-Parameter  impedi 

ance  model.  In  WSCAFFIP,  the  complex  ground  impedance  is  handled  by  the  graphical  user 
ini6ri3.cc. 


Table  2.  Estimation  of  Q  and  a  for  various  ground  surfaces 
Sf  =  0.875 


given  n'  =  0.750  and 


Description  of  Surface 

Q 

of  css) 

Dry  snow,  new  fallen  0.1  m  over  about  0.4  m  old  snow 

0.850 

23 

Sugar  snow 

0.825 

48 

In  forest,  pine,  or  hemlock 

0.825 

48 

Grass:  rough  pasture,  airport,  public  buildings,  etc 

0.675 

330 

Roadside  dirt,  ill-defined,  small  rocks  up  to  0. 1  m  mesh 

0.575 

960 

Sandy  silt,  hard-packed  by  vehicles 

0.475 

3470 

Clean  limestone  chips,  thick  layer  (0.01  to  0.025  m  mesh) 

0.425 

6470 

Old  dirt  roadway,  fine  stones  (0.5  m  mesh)  intstices  filled 

0.400 

7500 

Earth,  exposed  and  rain-packed 

0.350 

17100 

Quarry  dust,  fine,  very  hard-packed 

0.300 

41700 

Asphalt,  sealed  by  dust  and  light  use 

0.250 

120000 
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3.  Built-in  Features 


3.1  Sound  Speed  Parser 

As  discussed  in  Section  2.2,  the  FFP  uses  constant  sound  speed  layers.  In  previous  versions  of  the 
propagation  model,  a  novice  running  the  model  could  induce  error  into  the  computation  by  having 
the  atmospheric  layer  thickness  too  large.  Also,  previous  versions  of  the  propagation  model  used 
the  entire  meteorological  profile  regardless  whether  it  was  required  to  calculate  the  correct 
solution.  Although  this  did  not  induce  an  error,  the  run  time  was  slower  since  the  run  time  de¬ 
pends  on  the  number  of  layers  in  the  profile.  Thus,  a  sound  speed  parser  was  developed  to  miti¬ 
gate  these  two  problems.  The  parser  runs  invisibly  to  the  normal  user  and  yet  provides  a  sound 
speed  profile  with  the  appropriate  number  of  layers  to  the  maximum  height  required  to  give  a 
valid  solution. 

As  Figure  5  shows,  the  deviation  of  a  constant  sound  speed  layer  from  the  actual  value  of  the 
sound  speed  is  related  to  the  slope  of  the  sound  speed  profile.  The  larger  the  slope  of  the  profile, 
the  larger  the  difference  of  the  true  sound  speed  from  the  constant  in  the  layer.  We  can  minimize 
this  by  decreasing  the  thickness  of  the  constant  sound  speed  layer,  but  this  can  result  in  too  many 
constant  sound  speed  layers  in  the  regions  where  the  slope  is  very  small.  The  layer  thickness  also 
depends  on  the  wavelength.  The  shorter  the  wavelength,  the  smaller  the  layer  thickness.  By  using 
the  Green’s  Function  Parabolic  Equation  (GFPE)  model  [9]  as  a  comparison  base,  we  ran 


N_ 

Upward 

refraction 

1 

bJQ 

* 

<D 
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. .  .  _  1-  '  .  .  .  .. 

Sound  speed 

Figure  5.  The  difference  between  a  real  sound  speed  profile 
(solid)  and  WSCAFFIP’s  interpretation  (dashed) 
of  the  sound  speed  profile. 
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N^~ 

sound  propagating  above  the  duct  will  not  returnTn^therlg^01  ^  the  duCt’  any 

refracting,  This  means  that  usinpfhP^nf  r,  •  cethl ere§10n  above  the  duct  is  all  upward 

above  the  duct  and  not  affect  the  final  solufi™  c*  ^  W^Stefu1’  any  sound  above  tbe  duct  will  stay 
change,  with  azimuth  Z  ^  ^  ^ 

The  Simple  assumption"  ZZl  T  T  l ^  *°  so,™S «■*  Problem, 

targe,  and  Z  ^  ***  ^  S°“"d  «<>•«  the 

acoustic  energy  will  contr  ibut  t  th  P!i0n  ^  ^3Sed  °n  ray  pace  theory  that  suggests  that  no 
and  higher  is  lays  12 ta ^  Tt  ***  *  *  ^ 

fion  normally  works  very  well  but  if  thp  .  ,  .  eceiver  beight  of  interest.  This  assump- 

speed,  yet  still  smaller al^s  t  “  I  ^  18  C,°Se  t0  the  lvalue  of  the  sound 
contribute  to  the  solution  This  is  amanf  [T6  f? T ^  heightS  higher than  the  cutoff does 
theory  of  acoustics At  this  tin,  Zl  "  the  ray  and  -ve 

does  not  appear  very  often  For  insight  intoth  ^  I*1?11  ^  ^  probIem’  However,  it  apparently 
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Figure  7.  Sound  speed  profile  to  the  north,  east,  and  south  of  the 
target. 


3.2  Output  Smoothing 

The  attenuation  of  sound  with  range  appears,  as  a  simple  curve  with  range,  or  a  very  complex 
structure,  depending  on  the  number  of  paths  in  the  atmosphere  for  the  sound  to  reach  the  receiver 
and  on  the  various  modes  of  propagation  present  (surface  wave,  creeping  wave,  diffraction,  etc.). 
Figure  8  shows  an  example  of  the  complex  structure  of  the  attenuation  of  sound  with  range. 

When  experiments  are  typically  performed  in  the  atmosphere,  it  is  very  difficult  to  measure  the 
exact,  complex  structure.  This  is  because  the  atmosphere  is  dynamic.  Atmospheric  turbulence  is 
constantly  varying  the  temperature,  wind  speed,  and  wind  direction.  This  results  in  continually 
changing  fine  structure  in  the  attenuation  curve.  However,  since  acoustic  sensor  systems  record  a 
time  average  of  the  signal,  much  of  the  changing  fine  structure  will  average.  This  is  the  reason  for 
adding  the  output  smoothing.  The  model  performs  an  1 1 -point  average  on  the  attenuation  of 
sound  with  range.  This  acts  as  a  low-pass  filter  to  average  the  high  frequency  variations  in  the 
output.  The  high  frequencies  refer  to  the  oscillations  in  the  attenuation  curve  not  in  the  acoustic 
signal.  This  technique  provides  for  a  much  cleaner  graph  of  the  results  of  the  calculation.  Figure  9 
shows  the  same  case  as  Figure  8  except  with  smoothing  turned  on.  There  is  an  option  available  to 
switch  the  smoothing  off  if  one  is  interested  in  noting  the  difference  it  makes  in  the  quality  of  the 
output  graph.  Future  versions  of  our  models  that  will  include  turbulence  will  be  able  to  provide 
additional  information  about  the  bounds  of  the  amplitude  fluctuations  and  the  bearing  errors. 
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Figure  8.  Attenuation  versus  range  with  smoothing  turned  off. 
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Figure  9.  Attenuation  versus  range  with  smoothing  turned  on. 


3.3  Meteorological  Profile  Builder 

Since  the  atmosphere  plays  a  significant  role  in  affecting  the  propagation  of  sound,  knowing  the 
state  of  the  atmosphere  is  critical  in  performing  a  model  calculation.  Since  users  may  not  have 
meteorological  profiles  available  to  them  for  eveiy  location,  a  meteorological  profile  builder  was 
developed  to  give  the  user  a  tool  to  create  his/her  own  profiles,  based  on  basic  information.  The 
profile  builder  is  not  sufficient  to  generate  profiles  for  calculations  at  large  distances  (>  5  km). 
This  is  because  of  the  assumptions  upon  which  the  profile  builder  is  based.  For  long-range  propa¬ 
gation  calculations,  a  more  dynamic  meteorological  profile  model  must  be  used. 
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In  the  surface  layer,  mechanical  and  thermal  forcing  are  the  primary  driving  forces.  Monin 
Obukhov  introduced  scaling  parameters  that  allow  the  state  of  the  atmosphere  to  be  described, 
which  were  independent  of  height.  The  scaling  parameters  are  the  friction  velocity  (w*),  scaling 
temperature  (7*),  and  the  scaling  length  (X)  with  a  normalized  height  C,  =  z/L.  The  temperature 
and  wind  speed  profiles  are  calculated  by  (10] 


al¬ 

ii 

'a 

(  \ 

In 

z 

-Ym(z/L)+Wm(z0/L) 

K  \ 

T(z)  =  T(z,)  + 


7 IP, 


In  -  .^(z/D  +  TA/I) 
VV 


(3.1) 

(3.2) 


in  which  k  is  the  von  Karman’s  constant  (0.4),  z0  is  the  aerodynamic  roughness  length,  z,  is  an 
arbitrary  height,  and  P,  =  0.74.  The  universal  functions  used  in  this  development  are 


T,(0  =  n(0  +  V3 . 

<r<o 

(3.3) 

^m(C)=^(C)/Pt=- 4.7<r, 

<r<o 

(3.4) 

in  which 


«£)  = 


j(l-150~,/3 

C<  0 

11  +  4.7  C 

<r>o 

<r<o 

(3.5) 

[l  +  (4.7/ Pr)C 

o 

A 

The  problem  is,  what  are  good  values  for  the  friction  velocity,  scaling  temperature,  and  scaling 
length  for  different  locations,  time  of  day,  and  seasons  of  the  year?  In  order  to  determine  these 
parameters,  the  energy  balance  model  developed  by  Rachele  and  Tunick  [11]  was  employed.  The 
model  balances  the  energy  flux  at  the  soil  surfaces  as 

Rn=H  +  L'E  +  G  (3.6) 


In  which  RN  is  the  net  radiative  flux,  H  is  the  turbulent  sensible  heat  flux,  L'  E  is  the  latent  heat 
flux  attributable  to  surface  evaporation,  and  G  is  the  flux  of  heat  into  the  soil.  Each  term  in  the 
equation  reduces  to  the  basic  measurable  quantities.  RN  can  be  expressed  as  the  sum  of  the  fluxes 
of  the  net  shortwave  and  net  longwave  radiation.  The  energy  balance  model  is  constrained  to 
require  a  minimum  number  of  conventional  meteorological  inputs  from  a  2-m  reference  level.  The 
input  includes  temperature,  pressure,  relative  humidity,  and  wind  speed.  The  model  also  requires  a 
judgment  of  the  soil  type  and  moisture,  cloud  characteristics,  day  of  the  year,  time  of  day,  and 
longitude  and  latitude  of  the  site  of  interest.  Nearly  all  this  input  is  measured  for  most  areas  of  the 
world,  and  charts  are  available  to  estimate  the  soil  characteristics.  Appendix  C  provides  the  user 
with  a  series  of  tables  to  choose  some  of  the  parameters  that  cannot  be  directly  measured  such  as 
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albedo,  soil  characteristics,  etc.  Using  this  input,  the  energy  balance  model  calculates  the  similarity 

parameters  required  to  calculate  the  temperature  and  wind  speed  profile  required  by  the  acoustic 
propagation  model. 


3.4  Run  Options 

WSCAFFIP  offers  the  user  the  ability  to  set  program  flags  that  modify  how  certain  calculations 
are  made  and  generate  additional  output  files.  In  the  current  version,  there  are  several  run  options: 
source  to  receiver,  receiver  to  source,  output  sound  speed  profile,  disable  smoothing,  generate 
Battlefield  Acoustic  Simulation  System  (BASS)  output  file,  output  ray  information,  and  change 
the  units  of  some  of  the  meteorological  input.  The  default  configuration  is  for  source  to  receiver, 
output  sound  speed  profile,  atmospheric  pressure  in  millibars,  wind  direction  in  degrees,  and 
temperature  in  Celsius.  Here  is  a  description  of  each  run  option: 

source  to  receiver:  The  sound  is  attenuated  from  the  source  to  the  receiver.  This  is  used  when  the 
user  wants  to  view  the  sound  attenuation  from  a  source  (i.e.,  tank,  helicopter,  etc.) 

receiver  to  source:  The  sound  is  attenuated  from  the  receiver  to  the  source.  This  is  used  when  the 
user  wants  to  view  the  sound  attenuation  from  a  sensor  (i.e.,  WAM,  AHM,  etc.) 

output  sound  speed  profile:  Allows  the  user  to  output  the  sound  speed  profile  along  the  selected 
bearing  for  the  given  meteorological  profile.  This  option  generates  three  output  files:  sound.o, 
sound. new,  and  sound.inf.  The  first  file,  sound.o,  is  the  sound  speed  profile  calculated  from  the 
meteorological  profile  input  to  WSCAFFIP.  The  second  file,  sound.new,  is  the  sound  speed  profile 
used  by  WSCAFFIP  to  perfonn  its  calculations.  This  profile  was  generated  by  the  sound  speed 
parser  (see  Section  3.1).  The  last  file,  sound.inf,  was  also  generated  by  the  sound  speed  parser. 

The  file  gives  some  basic  analysis  of  the  sound  speed  profile  examined.  The  file  formats  for  the 
three  files  are  given  below: 

File  Name:  sound.o  (Original  Profile) 


Format: 

m 

C(l) 

Z(2) 

C(2) 

M 

M 

Z(Nint) 

C(Nint) 

in  which  Nint  is  the  number  of  interfaces,  Z(i)  is  the  height  of  the  ith  interface  in  meters  and  C(i)  is 
the  sound  speed  of  the  ith  interface  in  m/s. 
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File  Name:  sound. new  (Parsed  Profile) 

Format: 

Z(l)  C(l) 

Z(2)  C(2) 

M  M 
Z(N)  C(N) 

in  which  N  is  the  number  of  interfaces  in  the  new  profile. 

File  Name:  sound.inf  (Parsed  Profile  Analysis) 

Format: 

Results  of  Profile  Parsing 

Sound  Speed  at  the  surface:  342.42  m/s 

Sound  Speed  at  top  of  profile:  329.83  m/s 

Height  of  profile:  1327.00  m 

Largest  sound  speed  342.42  m/s  occurred  at  0.00  m 

Sound  speed  at  source:  340.58  m/s  Sound  speed  at  sensor:  342.42  m/s 

Highest  ray:  30.00  m 

Top  of  new  profile:  382.00  m 

Detector  just  outside  the  atmospheric  duct. 

Maximum  height  was  reset  to  capture  leakage  from  the  duct. 

Duct  Height:  382.00  m  with  sound  speed  of  338.22  m/s 
Possible  Sound  Ducts  Found 
Number  of  ducts:  3 

Duct  \#1  -  c:  341.955361  z:  172.000000  dz:  172.000000  index:  7 
Duct  \#2  -  c:  338.222784  z:  382.000000  dz:  182.000000  index:  13 
Duct  \#3  -  c:  338.158191  z:  450.000000  dz:  50.000000  index:  15 
Number  of  layers  in  original  profile:  27 
Number  of  layers  in  new  profile:  741 
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disable  smoothing:  As  discussed  in  Section  3.2,  WSCAFFIP  smoothes  the  higher  frequency 
variations  in  the  attenuation  table.  This  option  allows  the  user  to  switch  the  smoothing  off  if  the 
case  requires  it. 

generate  BASS  output  file:  If  the  user  needs  to  provide  attenuation  tables  for  the  older  BASS 
sensor  performance  model,  selecting  this  option  will  generate  an  additional  file  named 
proptab.dat.  This  file  has  the  format  of 


Format: 

dB(l,l) 

dB(2,l) 

L 

dB(M,l) 

dB(2,l) 

dB(2,2) 

L 

dB(M,2) 

M 

M 

0 

M 

dB(l,N) 

dB(2,N) 

L 

dB(M,  N) 

in  which  N  is  the  number  of  range  points,  M  is  the  number  of  frequencies,  and  dB(j,i)  is  the 
relative  sound  level  at  the  ith  range  and  jth  frequency  in  decibels  (re:  20  pPa). 

output  ray  information:  As  part  of  the  sound  speed  parser,  a  ray  trace  calculation  is  made  to 
determine  some  characteristics  of  the  atmosphere.  The  rays  generated  by  the  routine  can  be  sent 
to  a  file  named  “rays. out”  if  the  user  would  like  to  use  this  information. 

pressure  in  millibars:  The  current  version  has  changed  the  default  units  used  for  the  atmospheric 
pressure.  WSCAFFIP  now  defaults  to  use  millibars  as  the  measure  of  atmospheric  pressure.  To  be 
compatible  with  older  meteorological  profiles,  this  flag  was  provided.  If  the  user  is  using  an  older 
meteorological  profile,  which  employs  atmospheres  as  the  measure  of  atmospheric  pressure,  just 
unset  this  flag  to  correctly  read  the  file.  This  change  was  made  to  better  match  the  pressure  units 
normally  measured  by  meteorological  teams  in  the  battlefield. 

WSCAFFIP  uses  an  array  of  flags  to  internally  handle  the  run  options.  This  array  of  flags  can  be 
seen  in  the  configuration  file.  The  first  line  in  the  configuration  file  is  a  “printout”  of  the  flag  array. 
In  the  current  version,  there  is  allocation  for  as  many  as  16  run  flags;  however,  only  eight  flags  are 
currently  used.  The  remaining  eight  flags  are  for  future  use.  A  0  indicates  the  flag  is  off,  and  a  1 
indicates  the  flag  is  on.  The  order  of  the  flags  is 

Flag  Number  Function 

0  output  sound  speed  profile 

1  wind  direction  in  degrees 

2  source-to-receiver/receiver-to-source 

3  reserved 

4  reserved 
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5 

reserved 

6 

reserved 

7 

reserved 

8 

disable  smoothing 

9 

generate  BASS  output  file 

10 

pressure  in  millibars 

11 

temperature  in  Celcius 

12 

output  ray  information 

13 

output  ray  information 

14 

reserved 

15 

reserved 

4.  Operations 


4.1  Installing  WSCAFFIP 

There  is  little  to  do  for  the  installation.  WSCAFFIP  typically  is  sent  as  a  self-extracting  archive. 
The  user  runs  the  progam.  The  extractor  will  ask  for  the  directory  in  which  to  place  to  extracted 
files.  If  the  installation  is  a  regular  pkzip  file,  use  pkunzip  to  extract  the  files  into  the  directory 
where  the  program  will  operate.  The  archive  contains  four  files:  wscaffip.exe,  proptab.pth, 
proptab.ptd,  and  wscaffip.hlp,  with  a  sample  directory  containing  the  files  for  Appendix  A.  These 
files  will  be  extracted  to  the  directory  the  user  has  created  to  store  the  program.  The  proptab.pth 
and  proptab.ptd  files  contain  the  header  information  and  attenuation  table.  The  wscaffip.hlp  file 
contains  some  basic  help  information.  These  three  files  must  be  situated  within  the  WSCAFFIP 
directory  for  the  program  to  work  correctly. 

4.2  How  to  Run  WSCAFFIP 

When  the  user  runs  WSCAFFIP,  the  program  will  display  the  window  shown  in  Figure  10.  All  the 
user-required  information  is  given  under  the  “Model  Inputs”  menu.  There  are  five  selections  to  be 
made:  meteorological  profile,  geometry,  target  frequency,  ground  characteristics,  and  run  options. 
When  the  user  runs  the  program,  a  default  case  initializes  the  propagation  model.  The  user  can 
reset  the  program  back  to  the  run-time  defaults  by  selecting  the  default  settings.  Appendix  A 
shows  the  default  settings. 
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Scanning  Fart  field  Program  jSrAinP?  ; 


aeJcUnpUj  •  J?*  tjeuAc  Model  H^r-  • 
Ficquency 


0 


Range  [mj 


Figure  10.  WSCAFFIP  run  time  window. 


The  meteorological  profile  input  gives  the  user  three  options  from  which  to  choose:  open  profile, 
use  a  standard  meteorological  case,  and  create  a  meteorological  profile.  The  open  profile  is  used 
when  a  meteorological  profile  is  available.  The  format  of  the  meteorological  file  is  listed  below: 

Format: 

Nint 


Z(l) 

T(l) 

Rh(l) 

P(l) 

Wvel(l) 

Wdir(l) 

Z(2) 

T(2) 

Rh(2) 

P(2) 

Wvel(2) 

Wdir(2) 

M 

M 

M 

M 

M 

M 

Z(Nint) 

T(Nint) 

Rh(Nint) 

P(Nint) 

Wvel(Nint) 

Wdir(Nint) 

in  which  Nint  is  the  number  of  interfaces,  Z(i)  is  the  height  of  the  ith  interface  in  meters  also  Z(l) 
=  0  and  Z(i)  <  Z(i+1)  always,  T(i)  is  the  temperature  in  Kelvin,  Rh(i)  is  the  relative  humidity  in 
percent,  P(i)  is  the  atmospheric  pressure  in  millibars,  Wvel(i)  is  the  wind  speed  in  meters  per 
second,  and  Wdir(i)  is  the  wind  direction  in  degrees,  from  which  the  wind  is  blowing  relative  to 
north.  Each  meteorological  value  is  measured  at  height  Z(i).  Choosing  the  open  profile  selection 
displays  the  open  meteorological  profile  window,  as  shown  in  Figure  1 1 .  The  selected  profile  will 
be  read  and  used  to  perform  model  calculations. 
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The  next  way  of  obtaining  a  meteorological  profile  is  to  select  the  standard  meteorological  cases: 
overcast/calm,  overcast/windy,  sunny/calm,  sunny/windy,  night/calm,  and  night/windy  (see  Figure 
12).  The  six  cases  represent  typical  cases  that  the  user  would  observe.  Once  the  user  selects  a 
standard  case,  s/he  has  to  enter  the  temperature,  relative  humidity,  and  wind  direction  to  complete 
the  required  input.  If  the  user  would  like  to  further  customize  the  values  after  pressing  OK,  s/he 


r- Choose  one  case  - 


•v/V 


yF  -r.'v  Overcast,  Oaim  C*  Sunny,  Windy  i 


j ;Py  Overcast,  Windy  C"  Night,  Calm 


..jy ;  Sunny,  Calm 


,  :C:  ■ ...  Night,  Windy  \ 


T emperaturc  @  2m 
Relative  Humidity  @  2m 
Wind  Direction  @  2m 


+20.0 


25.0 


0.0 


degC 


Cancel 


Figure  12.  Standard  meteorological  profile  window. 
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can  select  “create  a  meteorological  profile.”  The  setting  for  the  standard  profile  case  will  be 
copied  to  the  “create  profile”  windows.  To  use  the  “create  profile”  windows,  refer  to  Section  4.3. 
Selecting  the  OK  button  results  in  the  creation  of  a  meteorological  profile. 

Next,  choose  the  target  frequency  input  selection.  There  are  three  choices  to  specify  the  frequency 
parameters  (see  Figure  13):  starting  frequency,  ending  frequency,  and  frequency  resolution.  The 
starting  frequency  is  the  first  frequency  used  in  the  frequency  sweep.  The  ending  frequency  is  the 
cut-off  frequency  for  the  sweep.  The  frequency  resolution  is  the  interval  between  frequencies  in 
the  sweep.  The  current  version  of  WSCAFFIP  does  not  allow  for  the  frequency  resolution  to  be  0. 
If  the  user  wants  to  run  one  frequency,  s/he  must  input  the  same  value  for  the  starting  and  ending 
frequencies  and  give  a  non-zero  positive  number  for  the  resolution. 


Figure  13.  Frequency  input  menu. 


The  next  input  selection  is  the  ground  characteristics.  The  options  for  this  are  shown  in  Figure  14. 
WSCAFFIP  gives  the  user  several  different  types  of  grounds  and  ground  conditions  from  which  to 
choose.  The  selection  here  will  generate  a  value  for  the  ground  impedance  used  by  the  propaga¬ 
tion  model;  see  Section  2.4  for  further  information.  A  future  version  will  allow  the  user  to  specify 
parameters  for  the  ground  impedance  model. 


to  Scanning  Fast  Field  Program  (SCAFFIPJ 


)  File  Configuration  I  Model  Inputs  .  Plot  Execute  Model  Help 
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Figure  14.  Ground  impedance  input  menu.  i 

The  last  input  selection  is  the  run  options.  The  options  for  this  are  shown  in  Figure  15.  These 
input  parameters  allow  the  user  to  modify  how  the  model  runs  and  generate  additional  output 
information.  The  details  of  the  options  and  the  default  values  are  discussed  in  Section  3.4: 


’ta  Scanning  Fast  Field  Program  (SCAFFIPJ 
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-  v  0 utputS bund  Speed  Profile 
Doable  Smoothing  /; 

*;:  £enei ate  BASS  Output  File 
.  Output  Ray  information./ 

v  Pressure  in  Millibar? 

•  v  Wind  Direction  in  Degrees'. 
V  Temperature  in  Celsius  • 


Figure  15.  Frequency  input  menu. 

Once  all  the  input  is  given  under  the  five  sections,  the  user  is  now  ready  to  run  the  propagation 
model  to  generate  an  attenuation  table.  To  do  this,  s/he  must  select  the  execute  model  menu  item. 
This  will  display  a  window  (see  Figure  16),  notifying  the  user  that  it  may  take  some  time  to 
perform  a  calculation  and  asking  the  user  whether  to  continue.  If  the  user  presses  the  “yes” 
button,  the  propagation  model  will  execute.  The  length  of  time  to  perform  a  calculation  varies, 
depending  on  the  height  of  the  meteorological  profile,  number  of  interfaces  in  the  profile,  number 
of  frequencies,  frequency  values,  and  propagation  range.  Run  tables  are  provided  with  each 
sample  case  in  Appendix  A.  The  current  version  of  WSCAFFIP  may  slow  or  lock  the  windows 
until  the  calculations  are  complete.  Future  versions  will  provide  an  indication  that  the  program  is 
still  executing. 
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Figure  16.  Model  execution  warning. 

4.3  How  to  Run  a  Meteorological  Profile  Builder 
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Figure  17.  Meteorological  profile  builder  window. 
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To  alter  these  assumptions,  the  user  must  select  the  advanced  button  on  the  “build  a  meteorologi¬ 
cal  profile”  window.  This  creates  another  window  shown  in  Figure  1 8.  This  window  contains 
additional  parameters  that  the  user  can  use  to  customize  the  meteorological  profile.  Most  of  the 
parameters  on  this  window  require  the  user  to  have  advanced  knowledge  about  the  environment. 
To  aid  the  user  in  selecting  values  for  this  window,  a  series  of  tables  has  been  provided  in  Appen¬ 
dix  B.  The  tables  will  allow  the  user  to  select  the  values  required  to  customize  the  profile  builder 
to  his/her  scenario.  Once  all  the  parameters  have  been  entered,  the  user  generates  the  profile  by 
selecting  the  OK  button  in  each  of  the  profile  builder  windows. 
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Figure  18.  Advanced  meteorological  profile  builder  window 


The  generated  profile  is  not  automatically  saved  to  a  file.  It  resides  in  memory  by  default.  It  is 
recommended  that  the  user  employ  the  save  option  under  the  model  input/meteorological  profile/ 
create  a  meteorological  profile  to  store  the  profile  to  a  file  for  later  use.  This  way,  the  user  will 
not  lose  the  profile  in  case  the  program  abnormally  terminates.  If  the  user  wishes  to  use  the  profile 
in  the  future,  s/he  just  has  to  read  the  profile  back  into  memory. 

4.4  Additional  Features 

In  addition  to  the  features  described  so  far,  WSCAFFIP  also  offers  the  user  some  additional 
capabilities:  store/retrieve/view  configurations,  save/retrieve  attenuation  table,  and  print  attenua¬ 
tion  plots  to  a  laser  printer.  The  program  has  been  tested  on  most  types  of  HP-compatible  laser 
and  inkjet  printers. 

In  order  to  preserve  the  user’s  settings,  WSCAFFIP  allows  the  user  to  save  the  current  input 
configuration  into  a  file.  Under  the  configuration  menu,  there  are  three  options:  open,  save,  and 
view.  Open  allows  the  user  to  open  and  read  a  previously  stored  configuration  file.  Save  will  store 
the  current  configuration  into  a  file.  The  view  option  allows  the  user  to  see  the  values  selected  for 
the  propagation  model  (see  Figure  19).  The  configuration  file  is  ASCII  with  the  format: 
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Geomeliy  - 
;  Source  Height  - 
f Receive*  Height 
[  Bearing. :  • 

Range.  -  J  ,  .  ’ 
i  •  R  arige  Resolution  : 


.1,00  m  . 
OibOdeg; 
P  40.00,00  m 
50.00(1(1  ' 


f ■  ground  Conditions  "  •  •  —  -  - 
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i  Omega.#!  ■ 
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Figure  19.  View  configuration  window. 

Format: 


I  I  Configuration  File  for  WSCAFFIP  v3.0 
//Run  Options 
options 
//  Geometry 

//  Source  Height  (m)  Receiver  Height  (m)  Bearing  (deg) 
Zsrc  Zdet  q 

"  Range  (m)  Range  Resolution(m) 

RDR 

//Frequency 

//  Start  Freq  (Hz)  Stop  Freq  (Hz)  Freq  Resolution  (Hz) 

AflVf 

//  Ground 

11  Number  of  Layers  Layer  Depth  (???) 

Ngnd  dgnd 

//Sigma]  Sfl  Pnl  Omega  1 
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s  1  sfl  Pn\  W1 

//  Sigma2  Sf2  Pn2  Omega2 

si  sfl  Pnl  W2 

//  Ground  Flag 

gndflag 

//  Meteorology 

//Meteorological  Profile  Options 

metflagl 

metflag2 

//  Size  of  File  Buffer 
buffer 

//Name  of  Meteorological  Profile 

metname 

//Plot  Title 

title 

in  which  run  options  are  the  1 6  binary  flags  defining  what  options  have  been  set.  See  Section  3.4 
for  a  list  of  the  options.  Zsrc  is  the  height  of  the  target  in  meters,  Zdel  is  the  height  of  the  detector 
in  meters,  q  is  the  bearing  in  degrees,  R  is  the  maximum  range  of  interest  in  meters,  D  R  is  the 
range  step  in  meters,/]  is  the  starting  frequency  in  hertz,/  is  the  ending  frequency  in  hertz,  D /is 
the  frequency  steps  in  hertz,  Ng#K/  is  the  number  of  ground  layers  (0  or  1),  dgn(j  is  the  depth  of  the 
top  layer  in  meters,  5']  is  the  flow  resistivity  in  cgs  units  of  the  top  layer,  sf  is  the  pore  shape 
factor  of  the  top  layer,  Pn\  is  the  grain  shape  factor  of  the  top  layer,  W ]  is  the  porosity  of  the  top 
layer,  $2  is  the  flow  resistivity  in  cgs  units  of  the  bottom  layer,  sf2  is  the  pore  shape  factor  of  the 
bottom  layer,  Pn2  is  the  grain  shape  factor  of  the  bottom  layer,  W2  is  the  porosity  of  the  bottom 
layer,  gndflag  is  the  type  of  ground  surface  used  in  calculation,  metflagl  is  the  meteorological 
profile  option  selected,  metflag2  is  the  flag  indicating  that  the  meteorological  profile  was  read 
from  a  file,  and  metname  is  the  name  of  the  meteorological  profile  used  for  the  calculation.  If  the 
profile  was  generated  by  the  profile  builder,  this  field  will  be  blank. 

4.5  Format  of  Output  Files 

Two  output  files  are  produced  by  the  program.  The  first  is  the  attenuation  table  called  proptab.ptd 
that  is  an  ASCII  table  with  a  format  of 
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Format. 


0  freq(l)  freq(2)  L  freq(M) 

R(l)  dB(],l)  dB(2,l)  L  dB(M,l) 

R(2)  dB(2,l)  dB(2,2)  L  dB(M,2) 

M  m  mom 

R(N)  dB(l,N)  dB(2,  N)  L  dB(M,N) 


“i"  ,he  and 

file  is  a  header  file  called  nrontah  ml,  J  h  ,f  que"cy  ln  declbeIs  <rt:  20  hPa).  The  other  output 

fonua,  of  the  header  file  is  ?  P  Pr°V,dm8  *bo* the  table.  The 


Format: 


//  Propagation  Table  Header 

// 

//  Created  by  WSCAFFIP 

// 


//  ID  String 

Generated  by  WSCAFFIP 
//  Prop  Data  Table  name 
proptab 

//  Prop  direction  re  wind 
propdir 

//  Number  of  ranges  (rows)  in  prop  table 
nranges 

//  Number  of  frequencies  (columns)  in  prop  table 
nfreqs 

//Atmospheric  Conditions  File 
none 

//  c(z)  data  table 
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none 


//  Number  of  heights  (rows  in  c(z)  table 
0 

//  Zmax  (0  if  unknown) 

0 

in  which  proptab  is  the  name  of  the  file  containing  the  propagation  table  with  the  extension  .ptd, 
propdir  is  the  propagation  direction  relative  to  the  wind  in  degrees,  nranges  is  the  number  of 
ranges,  and  nfreqs  is  the  number  of  frequencies.  The  remaining  four  fields  are  not  implemented  in 
the  current  version  of  WSCAFFIP.  The  default  names  for  the  attenuation  table  and  header  are 
proptab.ptd  and  proptab.pth,  respectively.  These  names  can  be  changed  with  the  file/save  table 
option.  The  extensions  need  to  remain  the  same  in  order  to  be  compatible  with  input  into  other 
sensor  performance  programs. 


5.  Future  Revisions 


This  is  the  last  planned  major  revision  of  WSCAFFIP.  Planned  future  revisions  will  be  to  eliminate 
any  “bugs”  found  and  to  optimize  the  program  to  work  with  Microsoft  Windows.  The  future 
model  development  will  focus  on  the  parabolic  equation  (PE)-based  propagation  models  such  as 
the  Finite  Element  PE  (FEPE),  Crank-Nicolson  PE  (CNPE),  and  the  Green’s  Function  PE 
(GFPE).  The  PEs  are  more  robust,  allowing  for  incorporation  of  additional  environmental  effects 
on  sound  propagation.  The  PEs  were  not  initially  used  because  of  several  difficulties  in  the  opera¬ 
tion  of  these  types  of  models.  The  FEPE  and  CNPE  are  very  computationally  intensive,  limiting 
the  frequency  that  the  user  could  practically  run.  The  GFPE  solves  these  problems,  but  some  of 
the  additional  input  parameters  were  difficult  for  the  non  expert  user  set.  The  incorrect  parameters 
would  result  in  incorrect  results,  which  a  non-expert  in  propagation  models  could  interpret  as 
errors.  However,  the  newer  version  of  the  GFPEs  is  greatly  improved,  and  many  of  the  input 
parameters  can  be  parameterized  and  automated,  based  on  the  typical  information  that  the  user 
has  available.  The  GFPE  will  allow  for  incorporation  of  scattering  from  atmospheric  turbulence 
and  irregular  terrain.  The  first  version  will  be  named  “SCAPE”  (Scanning  Parabolic  Equation). 

The  user  interface  and  capabilities  will  be  similar  to  WSCAFFIP.  The  main  difference  between  the 
two  is  that  SCAPE  will  allow  the  user  to  calculate  an  attenuation  table  for  several  heights  simulta¬ 
neously  instead  of  only  one  height  like  WSCAFFIP.  The  first  version  of  SCAPE  software  is 
available. 
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Appendix  A.  Sample  Case 


This  appendix  provides  a  sample  case  against  which  the  user  can  compare  his  or  her  version  of  the 
code.  This  aids  in  determining  whether  the  code  was  installed  correctly  and  training  for  the  new 
user.  The  meteorological  data  and  configuration  files  for  the  case  were  provided  with  the  program 
code.  The  sample  case  will  also  be  used  to  illustrate  some  of  the  environmental  effects  on  sound 
propagation.  Since  the  meteorological  data  file  will  be  provided  with  the  program,  tabulation  of 
the  input  data  will  not  appear  in  this  appendix  as  in  previous  user  manuals. 

For  this  case,  the  meteorological  profiles  of  temperature,  wind  speed,  and  wind  direction  are 
shown  in  Figures  A- 1  through  A-3,  respectively.  This  surface  wind  is  light  from  the  northeast. 
There  is  a  low-level  wind  jet  (local  maximum  in  the  wind  speed)  at  100  m  and  the  wind  direction 
is  from  the  northwest.  The  effect  of  this  low-level  wind  jet  will  be  obvious  when  the  sound  speed 
profile  is  examined.  As  can  be  seen  in  the  wind  speed  profile,  there  is  a  shift  in  the  wind  direction 
at  the  top  of  the  temperature  inversion  from  the  north  to  the  southeast.  Figure  A-4  shows  the 
sound  speed  profile  in  three  directions:  north,  east,  and  south.  This  represents  propagation  up¬ 
wind,  crosswind,  and  downwind. 


Figure  A -2.  Wind  speed  profile. 
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the  north  (0  the  sound  speed  profile  to 

poor  propagation  conditions.  This  condition  extends  th  '8-  AZ  1C!"m8  Upward  refraclion,  thus 
upper  half  of  the  wind  jet  (100  to  250  ml  cans  d  d  he  nnddle  ofthe  wind  jet  at  100  m.  The 

the  increasing  sound  speed  ZlZt  n  Z  ‘  ZT*  refraCti°n  reSi°n’  “  * 

away  for  the  ground,  causing  rap'd  “  "  ™  f  ^  ap 

acoustic  wave  reaches  1 00  m  above  ground  leve^  fAGL^Ah”  ZZ  W‘'h  Un“' 1,16 
will  be  refracted  back  toward  the  ground  f  ■  '  °Ve  this  height,  the  acoustic  wave 

range,  the  acoustic  ^  at  a  se, 

the  attenuation  field.  Looking  to  the  east  the  Zt  inn  ^  UP’ resuMmg  in  a  large  decrease  in 
upward.  This  means  that  the  attenuation  in  this  d'  ,  “  ^  "°W  downward  refractive  instead  of 

sound  speed  profile  to  the  south  also  contains  d  10”  If  "0t  “  great  aS  "  was  t0  the  north.  The 

Unlike  the  profile  to  the  east,  the  downward  "f™  °Ver  tte  flrst  1 00 

the  east.  The  impact  of  this  is  seen  in  the  attenuation  plotf  "  "0l  35  “Ck  as  « ia  “ 

north  ofthe  helicopter  (see  Figure  A T, w  The  att““ati<>"  * 

traced  wave  returns  ”  ““  '°SS  *he  re‘ 

to  the  upward  refractive  layer  in  the  sound  sveed  ”  n  inCre3Se  “ the  attenuatl0n  is  attributable 
with  A-6,  it  can  be  seen  that  ttelnlToZt ?  ^  A-5 

This  is  because  ofthe  absence  ofthe  .  f  mng6  1S  mUch  less  in  the  east  than  the  north. 

shows  a  large  downward  refractive  regkTextendkg^ ^OO3^  rtT  S°Und  Pr°file  t0  the  eaSt 
the  downward  refractive  laver  in  the  P5  ■  •  8  111  'n  ll°  atniosphere.  The  advantage  of 

graphs  eventually  show  similar ll  **  a,'enUa,i0n  b“«  *he  *» 

Finally,  Figure  A-7  shows  the  T-  *  ”iveS  “  ^re  A-5. 

there  is  about  the  same  amount  of  ^e^t *"  "  tl’a' 

m  me  east,  but  the  attenuation  is  larger  at  greater 
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Figure  A-5.  Attenuation  with  range  to  the  north  Figure  A-6.  Attenuation  with  range  to  the  east 

of  the  target.  of  the  target. 


Figure  A-7.  Attenuation  with  range  to  the  south  of 
the  target. 

distance  than  the  east.  This  is  partially  because  of  the  height  of  the  acoustic  duct  and  the  thickness 
of  the  region  between  the  downward  and  upward  refractive  regions.  The  duct  to  the  south  is  only 
1 00  m  thick,  that  is,  only  about  three  wavelengths  at  20  Hz.  Also,  the  thickness  of  the  transition 
region  is  very  small,  so  the  diffraction  of  the  refracted  wave  in  this  regions  leads  to  the  acoustic 
duct  being  “leaky”,  The  term  leaky  refers  to  energy  that  is  lost  from  the  acoustic  wave  traveling 
inside  a  duct,  thus  decreasing  the  energy  of  the  wave  as  it  travels  in  the  duct. 

Figures  A-8  through  A-10  compare  the  attenuation  of  sound  at  20,  100,  and  200  Hz  in  the  three 
directions.  It  can  be  seen  in  Figure  A-8  that  the  acoustic  wave  initially  loses  the  most  energy  but 
recovers  it  when  the  refracted  wave  returns  from  the  downward  refractive  region.  Also,  the  sound 
speed  profile  to  the  north  contains  the  largest  downward  refractive  region  that  traps  the  most 
acoustic  energy  within  the  duct.  Figures  A-9  and  A-10  also  illustrate  these  points  for  the  100-  and 
200-Hz  signals.  These  frequencies  show  the  effects  of  the  upward  refractive  region  better  than  the 
20-Hz  signal. 


Appendix  B.  Meteorological  Profile  Builder  Look-up  Tables 


Table  B-l.  IR  emissivity 


Ground  Cover 

E 

Fresh  snow 

0.99 

Old  snow 

0.82 

Dry  sand 

0.95 

Wet  sand 

0.98 

Dry  peat 

0.97 

Wet  peat 

0.98 

Soils 

0.9  to  0.98 

Asphalt 

0.95 

Concrete 

0.71  to  0.9 

Tar  and  gravel 

0.92 

Limestone  gravel 

0.92 

Light  sandstone  rock 

0.98 

Desert 

0.84  to  0.91 

Grass  lawn 

0.97 

Grass 

0.90  to  0.95 

Deciduous  forests 

0.95 

Coniferous  forests 

0.97 

Urban  area  (range) 

0.85  to  0.95 

Table  B-2.  Surface  albedo 


Surface 

a 

Surface 

a 

Snow  and  Ice 

Urban  (people  influenced) 

Snow,  fresh  fallen 

75  to  95 

Road,  black  top 

14 

Snow,  thawing 

20  to  65 

Road,  stone 

15 

Snow,  old 

40  to  70 

Road,  dirt,  wet 

18 

Snow,  icy 

75 

Road,  dirt,  dry 

35 

Ice,  grey 

60 

Road,  clay,  wet 

20 

Ice,  white 

75 

Road,  clay,  dry  & 

Ice,  water  covered 

26 

Road,  asphalt,  wet 

10 

Ice,  light  snow  cover 

31 

Road,  asphalt,  dry 

15 

Ice,  porous  and  melting 

41 

Parking  lot,  black  top 

8 

Desert  shrublands,  snow  covered 

18  to  19 

Concrete,  new,  white 

37 

Conifer  forest,  snow  covered 

59  to  67 

Buildings 

9 

Mixed  forest,  50  cm  snow  covered 

20 

Developed  urban  area  (average) 

15 

Grasslands,  snow  covered 

46  to  50 

Roof,  thatched,  new 

20 

Crops,  snow  covered 

18  to  19 

Roof,  thatched,  old 

15 
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Tundra,  snow  covered 

59  to  67 

Roof,  tiled,  dirty 

8 

Soil  and  Rocks 

Crops,  Natural  Terrain,  and  Vegetation 

Soil,  dark,  wet 

8 

Fallow  field,  wet 

5  to  7 

Soil,  dark,  dry 

13 

Fallow  field,  dry 

8  to  12 

Soil,  light,  wet 

Spring  wheat 

10  to  25 

Soil,  light,  dry 

Winter  wheat 

16  to  23 

Dark  organic  soils 

Rice  paddy 

12 

Dark  gray  silt 

- - - - 

12 

Winter  rye 

18  to  23 

Red  soils 

17 

Forest,  coniferous 

Clay,  wet 

16 

Forest,  deciduous 

10  to  20 

Clay,  dry 

23 

Grass,  green 

26 

Clay  loam,  wet 

19 

Meadows,  green 

10  to  20 

Clay  loam,  dry 

13 

Coniferous  trees,  dormant 

12 

Lime 

45 

Deciduous  trees,  dormant 

12 

Gypsum 

55 

Tall  grass,  dormant 

13 

Sand,  white,  wet 

25 

Mowed  grass,  dormant 

19 

Sand,  white,  dry 

35 

Tundra 

15  to  20 

Table  B-3.  Thermal  properties  of  soils 


Material 

State 

Soil  Thermal 
Conductivity 
(x  Wm~x  x  K~l) 

Soil  Thermal  Diffusivity 
(x  m 2  x  hr1) 

Sand 

Dry 

0.277 

0.000900 

Sand 

2.0 

0.002520 

Sand 

2.2 

0.002120 

Dry 

0.264 

0.000720 

8.3%  moist 

0.586 

0.001188 

Dry 

0.000450 

Yolo  silt  loam 

0.44 

0.001860 

Yolo  silt  loam 

0.848 

0.001730 

Sandy  loam 

Dry 

0.23 

0.000576 

Sandy  soil 

Saturated 

2.2 

0.002660 

Sandy  clay 

15%  moist 

0.925 

0.001330 

Dry 

0.3 

0.000549 

Moist 

1.47 

0.003960 

Clay  soil 

Saturated 

1.58 

0.001840 

Clay  pasture 

— 

2.76 

0.004320 

Calcareous  soil 

43%  moist 

0.712 

0.000648 

Soil  (generic) 

Very  dry 

0.167  to  0.345 

0.00072  to  0.0011 

Soil  (generic) 

Wet 

1.26  to  3.35 

0.000144  to  0.00036 

Mud 

— 

0.000792 

Peat  soil 

Dry 

0.13 

0.000540 

Peat  soil 

Saturated 

0.85 

0.000432 

37 


Table  B-4.  Roughness  length  for  different  surfaces 


Type  of  Surface 

z»(m) 

Ice/Smooth  mudflats 

0.00001 

Dry  lake  bed 

0.00003 

Tundra,  snow  covered 

0.0001 

Tundra,  patchy  snow 

0.0003 

Calm  open  sea 

0.0001 

Desert,  smooth 

0.0003 

Farmland,  snow  covered 

0.002 

Bare  soil,  tilled 

0.002  to  0.006 

Prairie 

0.01 

Grass  and  trees,  mixed 

0.035 

Brush,  scrub  growth,  open 

0.16 

Brush,  scrub  growth,  dense 

0.25 

Forested  plateau,  level  topography 

0.7  to  1 .2 

Forested  plateau,  rolling  hills 

1.25 

Subtropical  savanna,  grass,  few  trees 

0.31  to  0.41 

Subtropical  savanna,  grass  w/many  shrubs 

0.51  to  0.61 

Coniferous  forest 

1.1 

Residential,  low  density 

1.1 

City  park 

1.3 

Urban,  buildings,  business  district 

1.75  to  3.20 
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